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ABSTRACT 

We investigate the properties of massive, dense clouds formed in a barred galaxy 
and their possible relation to star formation, performing a two-dimensional hydro- 
dynamical simulation with the gravitational potential obtained from the 2Mass data 
from the barred spiral galaxy, M83. Since the environment for cloud formation and 
evolution in the bar region is expected to be different from that in the spiral arm 
region, barred galaxies are a good target to study the environmental effects on cloud 
formation and the subsequent star formation. Our simulation uses for an initial 80 Myr 
an isothermal flow of non-self gravitating gas in the barred potential, then including 
radiative cooling, heating and self-gravitation of the gas for the next 40 Myr, during 
which dense clumps are formed. We identify many cold, dense gas clumps for which 
the mass is more than lO^Af© (a value corresponding to the molecular clouds) and 
study the physical properties of these clumps. The relation of the velocity dispersion of 
the identified clump's internal motion with the clump size is similar to that observed 
in the molecular clouds of our Galaxy. We find that the virial parameters for clumps 
in the bar region are larger than that in the arm region. From our numerical results, 
we estimate star formation in the bar and arm region by applying the simple model of 
iKrumholz fc McKed (|2005l ). The mean relation between star formation rate and gas 
surface density agrees well with the observed Kennicutt-Schmidt relation. The SFE 
in the bar region is three times smaller than that in the spiral region. This trend is 
consistent with observations of barred galaxies. 

Key words: galaxies: ISM - galaxies: star formation - ISM: clouds - galaxies: structure 
- galaxies: kinematics and dynamics 



1 INTRODUCTION 

Star formation (SF) is one of the key processes govern- 
ing the evolution of galaxies. Many observations of nearby 
disk galaxies indicate an empirical relation between gas sur- 
face density (Ejas) and star formation rate surface density 
(Esfr)- T his relation is ca l led the Kennicut t- Schmidt (K- 
S) rela tion (|Schmidt 1 1| 19591 ). iKennicuttI l| 19981 ) ). lBigiel et all 
|2008l) show that T,sFR is well correlated with the hydrogen 
molecular surface density {T,h2) for local spiral galaxies. 

Many papers have been devoted to understand- 
ing the physical reason behind the K-S relation (e.g. 
iKrumholz fc McKed l|2005h ). However, it remains poorly un- 
derstood. Since massive stars are mainly observed in the 
spiral arms of disk galaxies, it has been proposed that star 
formation is regulated by galactic shocks driven by spiral 



density waves (e.g. iBinnev and Tramaine I (|2008|)) and b y 
the inc r ease o f cloud-cloud collisions jTasker fc 
[Tasked l|201ll )). IKrumholz fc McKed (|2005l ) propose a star 
formation model in which star formation efficiency depends 
on the turbulent properties of the molecular clouds and 
demonstrate that the Kennicutt-Schmidt relation can be ex- 
plained by their model. Turbulent internal motion in clouds 
can be exited by cloud-cloud collisions via the conversion of 
the orbital energy of th e clouds into i nterna l motion energy 



Bonnell et 



mer gy oi t h e clouds mto i nterna l motion energy 
all I2OO6I), iDobbs et all (|2006l ). iTasker fc Tainl 



20091) and [Tasked l|201ll )") 



Barred galaxies show different star formation activity 
in the bar regions and in the spiral arms regions, even 
when gas surface den si ty in both reg i ons is c omparable (e.g., 
iDownes et ahl (Il996l ). ISheth et al.1 l|2000l ). iMuraoka etall 
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l|2007ll , and iMomose et all (|2010l ')V Studying the physical 
reason for this discrepancy is an important step to under- 
standing the physics for the K-S relation. 

It is well known that the star formation activity is 
higher in spiral arms regions than in the bar regi ons, even if 
the bar regions are gas rich (Momoso et al. 2010). We focus 
on the possibility that this difference is related to a differ- 
ence in cloud properties between these regions, since the 
cloud environment is thought to be closely related with SF 
l|Krumholz fc McKed (^005)). 

Gas flows in the bar regions are elongated with strong 
shears and dark lanes th a t are evidence for strong shock 
waves( e.g. lWada fc Habd (| 19921 )). It is therefore natural to 
expect clouds forming in this environment to have different 
properties from elsewhere in the disk. 

In this paper, we present a numerical simulation of gas 
flow at high resolution in a barred galaxy potential. We in- 
clude cooling and heating processes and study the properties 
of clouds in different galactic environments within the disk 
and how this relates to star formation. For the purpose, we 
resolve gas clumps as large as molecular clouds. 

As the first step in our study, we perform a two 
dimensional simulation with spacial resolution of 4 pc. 
The galaxy model is that of a barred galaxy with 
a bar similar to M83. M83 is one of the best tar- 
gets to study the spacial variation of star formation 
efficiency in a barred galaxy, since it is nearby and 
has been well observed in various wave le ngths, e.g., 
atomic ga s llHuchtmeier fc B ohncnstGngcl''l98ll ), the molec - 
ular gas (iLundgren et al.l (12004a'). Sakamoto et al.l ll2004h 
Muraoka et al. | (|2007)), optical e mission lines ( Dopita et al.l 



2010l ). X-rav (Soria fc Wu]|2003l ) 



From the physical properties of the identified clumps, 
we estimate the SFE and SFR, using the simple 
star formation model of turbulent clouds proposed by 
iKrumholz fc McKed (|2005h . 

The plan of this paper is as follows: Section 2 describes 
our model and numerical method. In section 3, we show our 
numerical results. In section 4, we estimate SFR and SFE in 
th e bar region and the arm reg ion by using the simple model 
of IKrumholz fc McKed ([2005). In section 5, we present our 
discussion and conclusion. 



2 MODEL AND NUMERICAL METHOD 
2.1 Model Galaxy 

We use a gravitational potential for a barred galax y simi- 
lar to that of M83. This model is from the work of iHirotal 



20091) who analyz ed the 2Mass K-band image of M83 



Jarrett et aLll2003l ), assuming a constant mass to light ra- 



tio. The model galaxy consists of stellar bulge, stellar bar, 
stellar disk and dark halo components. The circular rotation 
velocity of the gravitational potential that is shown in Fig- 
ure [T] roughly_agrees_withtlM of molecular gas 
in M83 bv iLundgren et al.l (|2004bl ). The observed rotation 
velocity of molecular gas is 150 km/s at 2 k pc from the cen- 
ter an d 180 km/s at 3 kpc from the center (jLundgren et al.l 
l2004bl ). In our simulations, we assume that initial gas ro- 
tates in the gravitational potential with a circular velocity 
that balances the axial averaged gravity of the barred galaxy. 




Figure 1. The circular velocity in the model galaxy and the 
pattern speed of the bar. 



The pattern speed of the bar potential is Q,p — 54 km s ^ 
kpc~^, in keeping with the ob served global c haracteristics 
of the gas distribution in M83 (|HirotJ(|20o3), IffirotaeEaD 
(|2009l) ). The parameters for the initial conditions of our sim- 
ulations are shown in Table [T] For the initial radial distribu- 
tion of gas mass surface density in the disk, we assume the 
Gaussian central co mponent and the expon ential gas disk 
component given bv ILundgren et al] (|2004ah (Fig. [5]). The 
total gas mass within a radius of 7.3 kpc is 3.5 x 10^ M©. 



2.2 Numerical Method 

We simulate two dimensional ga s flow in the model galax y 
using our M-AUSMPW+ code (|Namekata fc Habd l201ll ). 
which is based on a n advection upstre am splitting scheme 
(|Kim fc KimI l2005bl) with the MLP5 (iKim fc KimI l2005al ) 
for calculating the higher oder numerical fluxes. The size 
of simulation region is 12.6 kpc x 12.6 kpc and covers the 
whole stellar disk of M83. The grid size in our simulation is 
3125^, with a cell size of 4 pc. 

For the first 80 Myrs, we calculate gas motion assum- 
ing an isothermal, non self-gravitating gas with T — lO'^K. 
At this stage, gas density distribution is steady. After this 
time, we take into account radiative cooling, heating and 
self- gravity of the gas for the next 40 Myr, which is long 
enough for the formation of the molecular clouds. Self- 
gravity of the gas is calcul ated by FFT of which a de- 
tailed description is given in lNamekata fc Habd (|201ll ). We 
do not include star formation processes or any stellar feed- 
back throughout our simulation, as we are concentrating on 
the properties of the clouds formed in the bar galaxy po- 
tential. We use a cooling function of gas w ith a solar metal- 
licity given by ISpaans fc NormarJ (|l997l ) . For the heating 
processes, we assume a uniform FUV radiation field and 
a uniform cosmic ray heating in the energy conservation 
equation in the numerical hydrodynamics. We use the far 
ultra-violet (FUV) heati ng rate Fpuv = 1-0 x lO'^^eGo 
dGerritsen and Ickdll997l). We assume e = 0.05 and Go = 
1.0 ( Huchtmeier fc Bohnenstengell Il98ll ). When radiative 
cooling is allowed, the gas is allowed to cool to 10 K. For 
calculation of the cooling rate, gas density is obtained by 
assuming the thickness of the gas disk is 70 pc. 
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Table 1. Simulation condition 



simulation grid points 


initial mass 


initial mass 


distribution 


I 3125 X 3125 


3.4 X W'^Mq 


Gaussian + 


exponential 



Mq is the solar mass. 



120Myr - 
Gaussian + exponential - 




1000 2000 3000 4000 5000 6000 7000 8000 9000 
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Figure 2. The azimuthally averaged gas surface density at t = 
120 Myr (solid line). Dotted line is the initial gas distribution. 



2.3 Clump finding method 

In order to study the properties of clumps in the bar and 
spiral region, we identify clumps that are dense and have a 
low temperature by the following method: (A more detailed 
description is given in iNa mekata fc Hab e (,2011i lV First, cells 
with a surface density higher than Sci and temperature lower 
than Tci are selected as a 'candidate' clump member. Next, 
if a neighbouring cell to the 'candidate' is also a candidate, 
these cells become members of the same group. We iterate 
this procedure for all candidate cells. If cell number of a 
group exceeds rid, the group is identified as a clump. 

Since the mean surface gas density is less than 
GOMq/pc^ except for where r < 600 pc, as shown in Fig- 
ure m we assume Ed = 60Mq /pc^ and rid ~ 12. These 
values of Ed and rid correspond to the lower Umit mass of 
clumps of m-d ~ W^Mq. 



the dense gas in the barred region shows dust-lane like fea- 
tures that extend away from the center and follow the in- 
ner edges of the spiral arms, similar to those observed in a 
barred galaxy. These dense gas features have a low temper- 
ature (T < 100 K ) and a re similar to the prev ious studies 
(|Wada fc NormanI (l200ll ). IWada fc Kodal (|200ll )). From our 
numerical result as shown in Figure |3] we call a rectangular 
region ( — 2 kpc < x < 2 kpc and —1 kpc< y <1 kpc) the 
bar region and a ring region between r = 2 kpc and r =4 
kpc the spiral region. 

There are many dense gas cells with E > QOMqpc~^ 
and T < 100 K. Total number of the dense gas cells is 
341,854 at t = 120 Myr. 

We show the probability distribution function (PDF) of 
surface density in Figure [S] This plot excludes the galactic 
central region of r < 500pc where there is a large concentra- 
tion of gas. The Figure shows that the PDF rapidly changes 
from t =80 Myr to 100 Myr and is almost steady in t ^ 100 
Myr in the range of E >10 Mq/pc^. T he PDF in this stage 
is app ro ximated by a log-n o rmal f orm ( IPadoan fc Nordlundl 
('200^, IWada fc NormanI l|200lf )) in the range of E > 
IOOMq/pc^. The following analysis focusses on the simu- 
lation results at t = 120Myr, since the time between t = SO 
and t = 120 Myr is long enough to have allowed more than 
one free-fall time for a typical giant molecular cloud (GMC), 
e.g. tff = y/3n/32Gp = 4.74(nj//100cm-3)-° '^Myr, and 
current estimates s uggest that GMCs liv e betw een 1 and 2 
free-fall times (e.g. iMcKee and Ostriker I (|2007l ). and refer- 
ences therin). 

We show the radial distribution of the azimuthally av- 
eraged gas surface density at t — 120 Myr in Figure [21 
There are gas concentrations in the galactic central region 
(r < 500pc) and near the bar ends (r ~ 3 kpc). These char- 
acteristic features agree with the observed radial profile of 
CO ( J = 1 - 0) and CO (J = 3 - 2) in M83 l|Muraoka et all 
[2003), except for the central 500pc region. 



3 RESULTS 

3.1 Disk gas evolution 

Figure [3] shows the gas surface density distribution at the 
final stage of the isothermal gas simulation at t = 80 Myr. 
In this Figure, rotation of the gas and the bar pattern speed 
are in the counter-clock direction. The major axis of the bar 
is along the x-axis. Two dominant gas spirals extend from 
the bar ends and dense gas distributes in the bar region. We 
continue the simulation with cooling and heating processes 
and self-gravity of the gas, untill t = 120 Myr. 

Figure |4] shows the gas surface density distribution at 
t — 120 Myr. Many irregular structures of dense gas have 
formed in the spiral and bar regions. In Figure |4l dense gas 
extends in the barred region more than in Figure [3] and 



3.2 Physical properties of the clumps 

We select clumps from our numerical results by the method 
described in section 2.3. These clumps are in the mass range 
of 10* < M < 10^ M0. The most massive clump is in r < 100 
pc from the galactic center. The total mass of the clumps 
is 1.7 xlO'"*Af0 in the gas disk. The mass of the clumps in 
the bar region is 1.32xlO^M0 and in the arm region is 0.38 
xIO^Mq. The total number of clumps in the gas disk is 
3122. The number of the clumps in the bar region is 495 
and in the arm region is 2617. 

We show a mass function of clumps in our numerical 
result at f = 120Myr in Figure [B] In this Figure, we plot the 
number of clumps with an equal spacing of d logjQ M = 0.05. 
From FigureEl dN/dM oc M'^'^^ in M > 10^ M©, since 
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X [kpc] 

Figure 3. The gas surface density map at t = 80 Myr. Until 
tliis stage, gas is isotiiermal. Tiie color bar in the right hand side 
shows the logarithmic scale of the surface density of gas in unit 
of Mq pc-2 
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Figure 4. The same as figure|3] but for gas surface density map at 
t = 120 Myr. Our numerical simulation includes radiative cooling 
after t = SOMyr. 



(1) 



dN _ 1 dN 
dM ^ M dlnM' 

This power index agrees well with t he observed values ^ 
-1.5 in M > IO^Mq in our Galaxy (jSolomon et al.lll987^ . 
We test different thresholds of surface density for the clump 
finding method, and confirm that the power index of the 
clump mass function is almost independent from the choice 
of the threshold value. 

We plot clump mass vs. their size in Figure [T] Plus 
symbols show clumps in the bar region. We use the color 
contour map to show clumps in the arm region which are 



too numerous to represent as points. The color bar shows 
number of clumps in the cell given by dividing the plot frame 
into 200 X 200. Since the clumps have an irregular shape, 
in this Figure, we estimate clump size / as: 



I 



pc-' 



where S is area of the clump. Large size clumps of I > 100 
pc are in the bar region and in the galactic central region 
(r < 500 pc). From Figure [T] we find the relation 



(2) 



that is shown in Figure [T] 

In Figure O we plot the one dimensional velocity dis- 
persion, o-„, of the internal velocity of the clumps vs. their 
sizes. (Ti, is given by: 




cl i,i 



(3) 



from our numerical result, where Ed i,j 
gas surface density and the velocity of gas in the cell 
belonging to this clump, respectively. Vd is velocity of the 
center of mass of this clump. Relation of the velocity disper- 
sions and sizes of clumps is approximated as 



0.8 (l/lpcy-^km/s 



(4) 



the arm region, which is shown by a dashed line in 
Figure [8] This relation is simila r to the r esults of I LarsonI 
and lSolomon et al. J 19871'). Lars on' ("1981') gives (t„ = 
Solomon et a l. (1987) gives (j„ — 



098 



1.1 {l/lpcf-^^ km/s and 

1.0 ±0.1 {l/lpcf-^^°-°^ km/s. There are clumps with larger 
velocity dispersions in the bar region than in the arm region. 
The@large velocity dispersion of clumps in the bar region 
may be due to their formation process; that is, clumps form 
from cooled gas in gas flow with strong shear motion in the 
bar region. 

We plot internal velocity dispersion vs. mass of clumps 
in Figure [9] The relation between the internal velocity dis- 
persion and the mass of clump is approximated as 



''km/s. 



(5) 



The virial parameter is a useful measure of the gravita- 
tional stability of the clumps. Virial parameter of a spherical 
clump is given by 

5o"^,/ . . 

where Mmc is its mass, (t„ is its one dimensional velocity 
dispersion of internal motion, and I is its size. We modify 



since our numerical calculation is two-dimensional. In Figure 
1101 we plot the virial parameter a„ir vs. mass of clumps. Sev- 
eral clumps in the bar region have larger virial parameters 
than in the spiral region. Number of clumps with a^ir > 1 is 
476 of the total 495 clumps in the bar and 2273 of the total 
2617 clumps in the arm region. These gravitationally un- 
bound clumps would be transient, since a^ir > 1. Such tran- 
sient clumps are found in the numerical simulation of giant 
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Figure 5. The PDF of surface gas density in unit of Mq/pc ^ . 
We plot numerical results at t = 80, 90, 100, 110, and 120 Myr. 




Figure 6. The clump mass function at t = 120 Myr. 



Figure 7. Clump mass vs. clump size for the clumps in our nu- 
merical results at t = 120 Myr. Plus symbols show clumps in the 
bar region and the color contour map shows clumps in the arm 
region. The color bar shows number of clumps in the cell given 
by dividing the plot flame into 200 X 200. 




log,Q (clump size) [pc] 



mole cular clouds formation by spiral shocks ijPobbs et al.l 
l2006l 'l. In iDobbs et all (|2006l '). they show that a large part 
of the formed clouds are gravitational unbound and tran- 
sient. (a!® Mass fraction of clumps with a^jir < 1 is 6.26% 
of the total mass of clumps in the disk. 0.14% of the total 
mass of clumps is in the bar region and 6.1% of the total 
mass of clumps is in the arm region. Fraction of clumps 
with a^ir < 1 is much smaller in the bar than that in the 
spiral region. It is notable that clumps in the bar region are 
transient and the fraction of gravitationally unstable clumps 
is small, since these results can explain low star formation 
efficiency in the bar region. 

We plot galactic radial distributions of velocity disper- 
sion of clumps in Figure [11] and virial parameters of clumps 
in Figure 1121 There are more clumps with large virial pa- 
rameters in r < 2.5 kpc than that in r > 2.5 kpc. There are 
many clumps with a < 1 in the spiral region. From these re- 
sults, we expect many star forming clumps in these regions. 



4 STAR FORMATION 

We estimate the star formation rate from our numerical re- 
sults. Since the cell size of our numerical simulation, 4 pc. 



Figure 8. The same as Figure [J] but for the distribution of ve- 
locity dispersion of internal motion vs. clump size for the clumps 
in our numerical res ults at t = 100 Myr. T he dashed line shows 
the relation given bv lSolomon et al 



is not enough to resolve the star formation process by hy- 
drodynamical simulation, we use the star formation model 
proposed by Kru mh olz fc McKee (2005) as sub-grid physics. 
In their model, they assume that star formation occurs in 
Jeans unstable parts of the giant molecular clouds and is 
regulated by the turbulent motions in the giant molecular 
clouds. In their model, they assume that turbulent veloci- 
ties in the giant molecular clouds is described by a power 
law of the size of the turbulent motion, a probability density 
function for the density is a log normal type and star forma- 
tion occurs in the Jeans unstable part of the giant molecular 
cloud. They obtain a star formation rate for a giant molec- 
ular cloud as: 

Mgmc 



SFR = SF Ra- 



ta 



where 



37r 
32Gp 



= i.5r 



VlOOOcm-3 



-1/2 



Myr. 



(7) 



(8) 
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Figure 9. The same as FigurejT] but for the distribution of veloc- 
ity dispersion of internal motion vs. clump mass for the clumps 
in our numerical results at t = 120 Myr. 




Iogio(clump mass) [M. 



Figure 12. The radial distribution of virial parameter of clumps 
in the galaxy at t = 120 Myr. 



and 



SFRfi = 0.014 



/ Qvir \ 

V 1.3 J 



f-V 



(9) 



where M — CTu/cs and Cs is the sound speed of clouds. We 
estimate star formation rate in a clump as 



SFR = 



/ SFRft/i 



y ^dyn 



(10) 



where ttici is the clump mass, /mc is the fraction of molecular 
gas of it, and t^yn is given by 



" Tvck, " ^-^Klki^VO (sOAfopc-^ 
Ec! is the mean surface density of each clump. 



Myr.(ll) 



Figure 10. The same as Figure!?] but for the distribution of the 
virial parameter of clumps with clump mass at t = 120 Myr. 




1500 2000 
Radius [pc] 



Figure 11. The radial distribution of velocity dispersion of 
clumps in the galaxy at t = 120 Myr. 



4.1 Relation of Star Formation Rate Surface 
Density and Gas Surface Density 

We obtain a SFR from our numerical results by assuming 
that the SFR for each clumpsQis given by equation H10|) 
with /mc = 1, that is, the clumps are assumed to be mainly 
composed of H2 molecules. We show the SFR surface density 
vs. surface density of mass of clumps over the average scale 
of 500 pc in Figure [T3] (the arms region) and in Figure [T?] 
(the bar region). These figures show that the SFR surface 
density is larger in the arms region than the bar region for 
the same clump mass surface density. Correlation shown in 
Fi gure 1 131 and Figur e 1141 agrees well with th e relation given 
bv lKennicut'^ (|l998l ) and lBigiel erall l|2008l ). 

4.2 Radial Distribution of Star Formation Rate 

We show the radial distribution of the mean star formation 
rate surface density in Figure [TS] In Figure 1151 we exclude 
the central region of r < 500 pc, since there is high mass 
concentration of clumps and the star formation rate is very 
high in this region. The star formation rate between 0.5 
kpc < r < 2.5 kpc is Esfr ~ IQ'"^'^ Mq /pc? /yr which is 
smaller than ~ 1O~'^'*M0 /pc? /yr in the region of 2.5 kpc < 
r <4 kpc. There is an enhancement of the star formation 
rate near r =3 kpc that is near the outside edge of the bar 
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Figure 13. The distribution of the SFR surface density aver- 
aged in the scale of 500 pc with the clump mass surface density 
averaged in the same scale at t = 120 Myr in the arms region. The 
color bar shows number of cells of numerical hydrodynamic simu- 
lation. The red solid line shows the Kennicutt - Schmidt relation. 
The blue dotted line shows 10 times the Kennicutt - Schmidt 
relation. The magenta dotted line shows 1/10 the Kennicutt - 
Schmidt relation. 



Kennicutt-SchmidJ 
Kennicutt-Schmidt- ' 
Kenricutt-Sclimidt-1 




Figure 15. The radial distribution of star formation rate ob- 
tained by using the Kulmholts-McKee star formation model in 
the model galaxy . We show the star formation rate in r ^ 600 
pc. 
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Figure 16. The radial distribution of star formation efficiency in 
r > eOOpc. 



Figure 14. The same as in figure [131 but for the bar regions. 



re gion. The star format ion rate agrees with the estimation 
of iMuraoka et all l|2007l ) within their uncertainty. 

We show the radial distribution of star formation ef- 
ficiency, i.e., star formation rate per unit clump mass, in 
Figure 1161 In Figure 1161 we also exclude the central re- 
gion of r < 500 pc. The star formation efficiency in the 
bar region is three times smaller tha n the arm region. Thi s 
trend is observ e d in b arred galaxies r iMuraoka etall l|2007l ). 
iMomose et aD (|201G| )) . Our numerical results show that 
many more clumps in the bar region have large virial param- 
eters with Q > 1 than in the arm region. This is the reason 
why the star formation efficiency is small in the bar region. 
In NGC 4303, SFE is larger than our numerical results with 
~ 10"*'^ yr"^ in bar region and ~ 10~*-^-10"'^-^ yr"^ in the 
arm region. This is because SFR is slightly higher in NGC 
4303 than the mean value of the KS relation (jMomose et al.l 
I2OIOI ). 



5 DISCUSSION AND CONCLUSION 

Our numerical results have shown that most of massive 
clumps in the bar region have large virial parameters, a 
and that mass fraction of bound clumps in the bar region 
is smaller than in the spiral region, although the mean sur- 
face gas density is larger in the bar region than in the arm 
region. Because of large shear gas flow in the bar region, 
most of clumps formed from gas in the bar region may have 
large internal gas motions as shown in Figure [TT] This is the 
reason why the mass fraction of bound clumps in the bar 
region is smaller than in the spiral region. 

The star formation rate estimation from our numerical 
results agrees well with the K-S relation as shown in sec- 
tion 4. In this a nalysis, we use the star formation model of 
iKrumholz fc Mc Kee (2005) in which they assume that grav- 
itationally unstable parts of the molecular clouds with small 
turbulent velocities form stars. In their model, they assume 
the scaling rela tion between tu rbulent velocities and turbu- 
lent sizes as in I LarsonI (|l983). In their model, SF rate is 
given by the gravitational unstable part mass divided by its 
free fall time. In our analysis, turbulent velocities of clump 
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sizes are assumed to be the velocity dispersions of internal 
motions of clumps in our numerical results. 

We have examined the radial distributions of star for- 
mation rate and star formation efficiency, since our estima- 
tion of star formation rate agrees well with the observed 
Kennicutt-Schmidt relation. Both radial d istributions are 
within the error bars of observation of M83 l|Muraoka et al.l 
in r > Ikpc and are smaller in the bar region than in 
the arm regions. The decrease of both quantities in the bar 
regi on is similar to the r ecent observations of barred galax- 
ies (|Momose et al.ll2010^ . Our numerical results show that 
most of massive clumps in the bar region have large virial 
parameters, a and that mass fraction of bound clumps in the 
bar region is smaller than in the spiral region, although the 
mean surface g as density is larg er in the bar region than in 
the arm region. iLundgren et al. I (2004bl ) shows that molecu- 
lar gas velo c ity di spersion is large in the bar region in M83. 
ISorai et al.l l|2012l ') shows that molecular g as in the bar in 
Maffei 2 has also large velocity dispersion. iLundgren et al.l 
l|2004bl ) estimated the Toomre unstable criteria of molecular 
gas and showed a large gas mass ratio to the critical mass 
of the Toomre instability in the bar region. However, there 
are few HIT regions in the bar region, except for the central 
star burst region, while HII regions correlate with the value 
of the ratio in the outside of bar region. This means that 
Toomre instability argument is not sufficient to discuss SF 
rate in the bar region. 

Our numerical results show that the turbulent star for- 
mation model can explain the property of star formation in 
barred galaxies. It is important to study observationally the 
difference in cloud properties in the bar and in the arm re- 
gion, and the relation between cloud property and SF activ- 
ity in various environments in galaxies using radio telescopes 
with high spatial resolution, e.g. ALMA, that will make it 
possible to resolve molecular clouds in extragalaxies. 

In our simulation, we do not consider feedback from 
star formation. If we consider destruction of clouds by 
star formation, the distribution of dense clumps would be 
changed and concentrate near the dust lanes more than in 
the present numerical result. In order to consider feedback 
process from star formation, 3D numerical simulations are 
needed, since the destruction of molecular clouds and ex- 
pansion of gas via heating of newly formed stars and super- 
nov ae can be simula t ed by three dimensional simulations ( 
e.g. iTasker fc BrvanI l|200d )). 
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